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ABSTRACT 

We present a physical interpretation for the locations of the sources of radio 
emission in IM Pegasi (IM Peg, HR 8703), the guide star for the NASA/Stanford 
relativity mission Gravity Probe B. This emission is seen in each of our 35 epochs 
of 8.4-GHz VLBI observations taken from 1997 to 2005. We found that the mean 
position of the radio emission is at or near the projected center of the primary to 
within about 27% of its radius, identifying this active star as the radio emitter. 
The positions of the radio brightness peaks are scattered across the disk of the 
primary and slightly beyond, preferentially along an axis with position angle, p. a. 
= —38° ±8°, which is closely aligned with the sky projections of the orbit normal 
■^j- ■ (p. a. = — 49.5°±8.6°) and the expected spin axis of the primary. Comparison with 



simulations suggests that brightness peaks are 3.6^q 4 times more likely to occur 
"^1" ■ (per unit surface area) near the pole regions of the primary (latitude, |A| > 70°) 

than near the equator (|A| < 20°), and to also occur close to the surface with 
. ~2/3 of them at altitudes not higher than 25% of the radius of the primary. 
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Introduction 



IM Pegasi (IM Peg; HR 8703; HD 216489; FK5 3829) is the radio-bright binary star 
which served as the guide star for the Gravity Probe B (GP-B) mission, the spaceborne 
relativity experiment developed by NASA and Stanford University to test two predictions of 
general relativity (GR). This paper is the sixth in a series of seven describing the program of 
very-long-baseline interferometry (VLBI) undertaken in support of GP-B. In the first paper 



in th e series we give an introduction to GP-B and to this series (Paper I, IShapiro et al. 



201 ll ). In the second and third papers we report on the structure and its changes of each 
of the three extragalac tic reference sources, 3C 454.3, B2250+194, and B2252+172 (Paper 
II. lRansom et al.ll201ll ). and on the degree of astrometric stability of the "co re" of 3C 454.3 
in two extragalactic celestial reference frames (Paper III, iBartel et al.l 1201 if ). In the fourth 



and fifth papers, we describe our astrometric analysis technique (Paper IV, iLebach et al. 



2011) and present our result on the proper motion, parallax, and orbit of IM Peg (Paper V, 
Ratner et al.ll201ll ). In this paper (Paper VI), we discuss the locations of the sources of radio 
emission in the IM Peg system and give a physical interpretat ion of the sources' origins. In 
the last paper of our series (Paper VII, iBietenholz et al.ll201ll ). we discuss the radio images 
of IM Peg and include a movie of this star's changes over the duration of our 8.5 yr observing 
program. 

IM Peg is a close binary with orbi tal period ~24.65 days and an essentia lly circular orbit 
wi th an eccent ricity of 0.006 ± 0.007 fjBerdyugina. Ilyin. Sz Tuominen 1999 ). It is classified 



bv lHalll f 



also ESA 



9761) as an RS CVn. The system is at a di stance of 96.4 ± 0.7 pc (Paper V; see 



19971 ) and has an inclinatio n 65° < i < 80° (IBerdyugina et al.lll999t ILebach et al. 



19991 ). The primary is a K2 III star (IBerdyugina et al.l 119991 ) which is magnetically active, 
showing bright emission features (e.g., Ca II H and K, Mg II H and K, C IV) that are 
presumably produced by high-temperature species in its chromosphere and transi tion region 
( iHuenemoerder. Ramsey, fc Buzasilll990l ; iDempsey et al.lll996l ; lOlah et al.lll998l ). In addi- 
tion, Doppler optical images of the photosphere of the primary exhibit la rge, relatively dim 



"spot " regions, covering collectively >15% of the visible stellar surface (IBerdyugina et al. 
20001 ). The sun-like secondary is ~60 times l ess luminous in the optical than the primary, 
and has also been detected spectroscopically (IMarsden et al.l 120051 ) . 



Radio emission from IM Peg was first detected by ISpangler. Owen, fc Hulsei (119771). 
Since then, IM Peg has been include d in two radio surveys of RS CVn systems (IMorris fc Mute! 
19881 ; iDrake. Simon, fc Linskyj|l989l ). but few details of its radio properties are published. Its 
flux density at centimeter wavelengths h as ranged between ~0.2 mJy and ~80 mJy (Paper I; 



Lebach et al. 



1999 



Boboltz et al.l 120031 ). and can be highly variable on sub-hour time scales 



( ILebach et al.lll999l ). VLBI observations of IM Peg were made in the early 1990's as part 
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of an astrometric pro gram to link the Hipparcos optical reference frame to the extragalactic 



19991 ). but no image of IM Peg from this program 



radio reference frame ( ILestrade et al.lll99 
was published. 

Models for the microwave radio emission of RS CVn binaries suggest three possible 
source regions for the emission within the system: (i) magnetic-loop structures attached to 



one s t ellar component, namel y the active subgiant or giant, in the binary (e.g., iMutel et al. 



19851 ; iFranciosini et al.l Il999f); (ii) a joint magnetospher e for the two components of the 
binary (e.g., Uchida fc Sakural 1983 ; Ransom et al. 2002 ); and (hi) the region between the 
two components ( ILestradelll996l ). Multi-epoch astrometric VLBI observations can potentially 
distinguish am o ng the se scenarios. In the related case of the close binary in the Algol system, 
Lestrade et al.l ( 119931 ) were able, with astrometric VLBI from four epochs, to identify the 
cooler K subgiant star, and not its B dwarf companion or the intermediary region, as the 
likely source of the radio emission. Until now, no such identification has been made for any 
RS CVn system. 

A total of 35 additional sessions of astrometric VLBI observations of IM Peg were con- 
ducted between 1997 and 2005 in support of GP-B. Consequently, IM Peg is now more 
extensively observed at centimeter wavelengths than any other binary. In §[2] below, we give 
the previously determined orbital parameters of IM Peg, and set the stage for a discussion 
about the source region of the radio emission in the binary system. In § J3j we give an 
overview of the VLBI observations and astrometric analysis p rocedures. In §|U we summa- 
rize the astrometric results for IM Peg presented in iPaper VI . focusing in particular on the 
apparent orbit of the radio emission region and the distribution of the residuals on the sky. 
We discuss our results in §[5] and give our conclusions in §|6j 



2. Previously Determined Orbital Parameters of IM Peg 

Optical spectroscopic and photometric observations provide accurate values for most of 
the basic physical properties and orbital elements of the IM Peg binary system. Table [T] 
summarizes these results. The orbits of the primary and secondary stars projected on the 
sky are particularly relevant to the problem of determining the location of the radio emission 
within the binary system. If the emission source is closely tied to either of the two stars, 
then it likely travels nearly the same projected orbital path as that star. However, if the 
emission source arises primarily in the interbinary region, the source could remain more 
nearly stationary near the center of mass of the binary. 



The rotation of at least the K2 III primary is synchronized with the star's orbit ( lOlah et al 
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Table 1. Properties and Previously Determined Orbital Parameters of IM Peg 



Parameter 


Value 


Reference 3, 


Trigonometric Parallax (mas) 
Distance (pc) 


10.33 ±0.76, 10.370 ±0.074 
96.8;tg;2, 96.4 ±0.7 


1,2 
1,2 


Stellar Properties 6 


Mass (Mq) 
Spectral Type 

(K) 
Radius (Rq) 
Radius (mas) d 


1.8 ±0.2 1.0 ±0.1 
K2 III G V? c 
4550 ± 50 5650 ± 200 c 
13.3 ±0.6 1.00±0.07 c 
0.64 ± 0.03 0.048 ± 0.004 c 


3,3 
4,3 
4,3 
4,3 
4,3 


Orbital Elements' 3 


asini (Rq) 
asini (mas) d 
P (days) 

i n 

e 

T coni (HJD) e 


16.70 ±0.02 30.34 ±0.03 
0.806 1.464 
24.64877 ± 0.00003 
65. ..80, >55 
0.0 (assumed) 
2450342.905 ± 0.004 


3,3 
3 

4,5 
4 
3 



a First reference is for the first entry, second reference, if present, is for the second 
entry. 

b Two entries for lines 3—9 correspond to the two stars of the binary system, with 
entries for the primary listed first. 

c The spectral type, effective temperature, and radius of the secondary are inferred 
from the flux ratios (at two wavelengths) of the two stellar components and the values 
for the radius and effective temperature of the primary under the assumption that the 
secondary is a main sequence star. 

d Computed for a system distance of 96.4 ± 0.7 pc. The uncertainty in the asini 
value in Rq units is not propagated into mas, since the uncertainty in the inclination 
is the dominant source of error in any spectroscopic determination of the semimajor 
axis. 

c Heliocentric time of conjunction with the K2 III primary behind the secondary. 



References. - 1 Hippa rcos Catalogue jESAl ll997T) ; 2. VLBI jPaper Vh ; 3. 
Marsden et al 11200511 ; 4. iBerdvugina et alj jl99Sh (e = 0.006 ± 0.007); 5. 
Lebach et al.1 Jl999f) . 
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19981 ). which is expected to be circ ular based on tidal theory (|Zahnlll9770. and indeed de- 



term ined to be very nearly circular (jOlah et al.lll998t iBerdyugina et al 



1999; Marsden et al. 



20051 ). Projecting a circular orbit with system inclination i > 55° (see Tabled]) on the sky 
yields a highly eccentric, elliptically-shaped orbit. The semimajor axis lengths of the pro- 
jected (elliptical) orbits of the primary and secondary, a\ and respectively, are equal to 
those of the true (circular) orbits, which are constrained by spectroscopic data. Using the 
values from such data for a-\ sinz and a?smi in Tab le [IJ the narrow range of allowed incli- 
nation val ues 65° < i < 80° ( iBerdyugina et al.lll999l ). and a system distance of 96.4 ± 0.7 pc 
( IPaper VI ). we deduce that the semimajor axes of the orbits for the primary and secondary 
in angular units are a\ = 0.85 ±0.04 mas and a<i = 1.55 ±0.06 mas. These values imply that 
the maximum dimension of the projected orbit of each star is large enough to be detected 
with astrometric VLBI and that the projected orbits of the primary and secondary stars 
are clearly distinguishable. If the radio emission from IM Peg is indeed spatially associated 
largely or entirely with one of the two stars of the binary, then a projected orbit derived 
from our astrometric observations can distinguish between these two possibilities. 



3. Overview of Observations and Astrometric Analysis Procedures 

3.1. VLBI Observations 



Our 35 sessions of 8.4 GHz (A = 3.6 cm) VLBI observations each used a global VLBI 
array of between 12 and 16 telescopes. For each session, we interleaved observations of 
IM Peg with either two or three extragalactic referen ce sources nearby on the sky, so that we 



could employ the phase-referencing technique (e.g. , IShapiro et al.l 11979c iBartel et al.l Il986 



Lestrade et al.lll990l ; iBeasley fc Conwayill995l ) and determine an accurate astrometric posi- 
tion for IM Peg. For a full descri ption of our array, typical observing schedule, and data 
recording parameters, see lPaper III . Other aspects of our observing st rategy, e. g., s easonal and 
orbit-phase distributions of our observation sessions, are discussed in lPaper Vl and lPaper VII . 
Since the locations of the sources of radio emission in IM Peg are of special interest in this 
paper, we emphasize here that we took care in our scheduling to achieve an approximately 
even distribution of orbital phase, without strong correlations between phase and either year 
or season. 
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3.2. The Astrometric Technique 

The phase-referencing process by which we estimated the position of IM Peg for each 
observing session included both phase-delay fitting with a Kalman-filter estimator and phase- 
reference mapping (see Paper IV for details). This process improves our astrometric accuracy 
by allowing us to effectively model the contributions of the troposphere, ionosphere, and the 
extended structure of IM Peg's radio emission regions (see § 13.3 j) . in spite of the often low 
flux density (as low as ~0.2 mJy) of the stellar radio emission. The final image of IM Peg 
produced for each session is referenced to the quasi-stationary component CI, the 8.4 GHz 
"core," of 3C 454.3. The identification and stationarity of this reference point are discussed 
in Papers II and III, respectively. 



3.3. Choosing a Position for IM Peg at Each Epoch 



As mentioned in iPaper VI . our phase-referenced images reveal three general categories 
for the radio source structure of IM Peg: (1) single-peaked with peak located near the 
center of a marginally extended source; (2) single-peaked with peak located off center of 
an elongated source; and (3) double-peaked (or in one case apparently triple-peaked) with 
maximum separatio n between peaks of ~1.5 mas. An example from each of these categories 



is given m 



Paper V 



For the definition of the "observed" position for IM Peg we considered three possible 
choices: (i) the position of the brightness peak, interpolated between pixels, of the image for 
each epoch, (ii) the central position of the single elliptical Gaussian fit to the image for each 
epoch, and (iii) the central position of the single fit Gaussian for each single-peaked epoch 
and the position of the unweighted midpoint between interpolated brightness peaks for each 
multiple-peaked epoch. 

We then fit the astrometric model described in § 13.41 below to each set of positions. We 
obtained the best fit for choice (iii): the chi-square per degree of freedom for the resulting 
weighted least-squares fit was 30% lower than for the worst fit, that for choice (i). We 
therefore adopted set (iii) as our standard set of 35 IM Peg positions for the estimate of 
position at epoch, proper motion, parallax, and orbital parameters of IM Peg. Set (iii) 
represents in effect smoothed data, taking into account the extended structure of the source. 
For the study of the locations of the sources of radio emission in this paper, however, the 
distribution of the locations of emission peaks is more important, and therefore set (i) is 
used for such an investigation, as described below in § 14.31 
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3.4. The Astrometric Model 



We fit to the 35 positions of set (iii) a nine-parameter model (see also iPaper VI ) de- 
scribing the motion of the IM Peg radio emission on the sky. The model parameters are the 
five scalar parameters needed to specify IM Peg's position at a reference epoch, its proper 
motion, and its parallax, plus four more to specify the projection on the sky of its (as- 
sumed) zero-eccentricity orbit of known period. To maintain the linearity of the model in 
all unknown parameters, we let the orbit parameters for each coordinate (a and S) be the 
amplitudes of the sine (A sa , A s $) and cosine (A ca , A c $) terms in orbital phase. 



4. Results 
4.1. Astrometric Solution 



The full astrometric solution is presented in Table 3 of IPaper VI . In this paper we focus 
on the orbit. For the convenience of the reader and for easier comparison with previously 
determined orbital paramet ers, we reproduce in Table [2] the values for the four orbit parame- 



ters determined in IPaper VI . The estimated orbit and the 35 residual position determinations 
with their corresponding orbit-model-predicted positions are plotted in Figure [TJ The in- 
ferred direction of motion is counter-clockwise as indicated by the arrow. The size, shape, 
orientation, and "timing" of the orbit do not change significantly if we choose instead set (i) 
or set (ii) of astrometric positions described in § 13.31 The orbit is also robust against reason- 
able changes in either the set of epochs included in the fit or the addition to the astrometric 
model of a constant proper acceleration. A full discussion of our error analysis is given in 



Paper VI . The orbit of the secondary is obtained from the ratio of the component masses 
given in Table [TJ In Figure [2] we show an artist's three-dimensional rendition of IM Peg with 
the primary and secondary in their estimated orbits as seen from Earth. 

We emphasize three points from Table [21 each of importance for our study of the location 
of the sources of radio emission: 



1. Our estimate of the length of the semimajor axis of 0.89 ± 0.09 mas ([Paper VI ) is both 
statistically significant and consistent within la with that of the semimajor axis inferred for 
the IM Peg primary from optical spectroscopy (see §[2J Tabled]). 

2. The ratio of the length of the minor axis to that of the major axis of the projected orbit 
of 0.30 ± 0.13 is relatively small allowing for a relatively accurate estimate of the position 
angle (p. a.) of the ascending node, Q = 40.5° ± 8.6°, and hence the p. a. of the projected 
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Fig. 1. — The orbit (dashed ellipse) derived from the nine-parameter fit to the set of 35 
positions of the IM Peg radio source. The inferred direction of motion on the sky is counter- 
clockwise, as indicated by the arrow on the ellipse. The ascending node is given as an open 
circle on the orbit in the northeast. In other words, the part of the ellipse with the arrow 
on it is closest to Earth. The observed position for each epoch is plotted (with an open 
triangle) after subtraction of the estimated position at reference epoch, proper motion, and 
parallax. This position corresponds to the peak of a single Gaussian component fit to the 
source region or the midpoint between two (or three) local maxima in the source region. A 
solid line connects each observed position with the corresponding position indicated by a dot 
on the estimated orbit (see text). 
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Fig. 2. — Artist's three-dimensional rendering of the IM Peg binary system as seen from 
Earth. Grid lines are drawn every 0.5 mas. The primary is the larger red star with dark 
spots, while the secondary is the smaller yellow star. The projected orbit of the primary 
is the same as the inferred radio source orbit shown in Figure [U The size of the orbit of 
the secondary is computed from that of the radio orbit and the ratio of the component 
masses given in Table [TJ The diameters of the primary and secondary stars correspond to 
the nominal values given in Table [TJ The black dot at the center of the figure represents the 
center of mass of the binary system. The system is shown with the primary at its ascending 
node. The spin axis of the primary is indicated by the rod passing through the star's center. 
It is oriented parallel to the orbit normal (see §[5]). Note that the visible pole of the primary 
is near the south-east side of the stellar disk. 
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orbit normal, Q - 90° = -49.5° ± 8.6°. 

3. The time of conjunction, T con j = 2450342.56 ± 0.44 JD, implied by our orbit is also con- 
sistent within la with that derived for the primary from optical spectroscopy (see Tabled]). 



4.2. The Mean Location of the Source of the Radio Emission 



In lPaper VI we show that the residuals of our nine-parameter astrometric fit to set (iii) of 
the 35 positions of IM Peg scatter about a well determined orbit on the sky. How does this fit 
determine the mean location of the radio emission source with respect to the optical primary 
and secondary of the binary system? With a near zero-eccentricity orbit, the axial ratio of 
the projected radio orbit corresponds to a n inclination of 73° ± 8°. Combining our value for 
the inclination with the asinz estimate of iMarsden et al.l (120051 ) for each stellar component 
(see Table [1]) leads to semimajor axes of 0.84 ± 0.03 mas for the primary and 1.53 ± 0.06 mas 
for the secondary. Thus our estimated semimajor axis of the radio orbit of 0.89 ± 0.09 mas 
agrees well with that of the spectroscopic orbit of the primary, but differs significantly from 
that of the secondary. Furthermore, the time of conj unction, T mn ), o btaine d for the radio 



orbit is only 0.3 ± 0.4 d earlier than that found by IMarsden et al.l (120051 ) (see Table [T]). 
The estimated radio orbit is thus in the same phase within the error as the spectroscopic 
orbit. The above difference in T con j corresponds to a physical offset between the center of 
the primary and the mean position of the radio emission from our model orbit that is only 
0.12 ±0.15 times the radius of the primary and is not significantly different from zeroQ That 
is, the offset is with la accuracy likely to be distant from the center by less than 27% of 
the radius of the primary. We therefore have strong observational evidence to conclude that 
the active primary is the source of the radio emission in IM Peg. This result makes IM Peg 
only the second close binary system for which such an identi fication has unambigu ously been 
made, the other being the close binary in the Algol system (ILestrade et al.lll993l ). 



lr rhis result depends on the assumption that there is no offset between the center of the primary and the 
mean position of the radio emission that is constant or steadily increasing/decreasing over the 8.5 years of 
VLBI observations. Such an offset could of course not be determined in our fit since it would be absorbed 
in our position and proper-motion estimates. However, as we demonstrate in § 14. 3i the distribution of our 
position solutions well covers the disk of the primary. Therefore it appears that the emission locations are 
very closely linked to the primary and any constant or linearly changing offset of the cente r of this d istribution 
from the center of the primary is likely smaller than the radius of the primary (see also lPaper VI ). 
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4.3. Distribution of Position Residuals on the Sky 

We show in Figure [Sh the residuals of our set (iii) position solutions to our nine- 
parameter weighted least-squares astrometric fit. The residuals correspond to the positions 
plotted in Figure [1] after removal of the model orbit. In addition we plot the disk of the 
primary, placing its center at the origin of the diagram. There are two important features 
in the sky-distribution of the residuals: First, almost all of the residuals lie inside of the 
disk of the primary, with some residuals going only slightly beyond it. In fact, the 0.55 mas 
root-mean-square (rms) scatter (0.35 mas in a and 0.46 mas in 5) is almost equal to the 
angular radius of the primary of 0.64 ± 0.03 mas (Tabled]) but smaller than the semimajor 
axis of the orbit of the primary of 0.89 ± 0.09 mas (Table [2]) . With a relatively small mean 
standard error in the position determination of the stellar radio source of 0.07 mas in a and 
0.09 mas in 5, the scatter is largely intrinsic to the emission source^] Second, the residuals are 
scattered preferentially along a northwest-southeast oriented axis, approximately along the 
projected normal of the orbit. It is likely that this preferential elongation of the distribution 
of the residuals is also intrinsic to the emission source. The synthesized interferometer beam 
is by comparison more nearly oriented north-south (mean p. a. ~ —7°). Moreover, when we 
fit for the purposes of error analysis the same nine-parameter model to the positions of the 
secondary reference source B2250+194 (which is ~5 times farther away from 3C 454.3 than 
IM Peg), we found an approximately threefold smaller scatter in the rms residuals and a 
smaller correlation between coordinates. 

For the physical interpretation of the location of the source of radio emission it is more 
meaningful to study the distribution of the positions of the brightness peaks (set i) which 
indicate where the dominant part of the emission originates. We therefore added to the 
residual for each epoch the difference between the position of the interpolated brightness 
peak (from set i) and the position of the Gaussian or midpoint (from set iii) used in our 
astrometric fit. Use of these "modified" residuals (hereafter referred to simply as residuals) 
ensures that we are plotting the offsets of the interpolated brightness peak at each epoch 
from our estimated orbit (see Figure Et) . The rms scatter about the mean is 0.62±0.03 mas, 
and, as expected, somewhat larger than the rms scatter of 0.55 mas in Figure [3J2. 

We determined the 0.03 mas standard error for the above sample estimate of rms scatter 



2 The mean standard error in the position determinations is the root-sum-square of the mean statistical 
standard error associated with determining the position of the reference point in t he image of IM Peg at each 
epoch and an estimated ~0.06 mas astrometric error in each IM Peg coordinate (jPaper VI ). The astrometric 
error inclu des the rms 'jitter' of 3C 454.3 component CI relative to our extragalactic reference frame (see 



Paper 111 ) . 
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Fig. 3. — (a) Sky distribution of the observed positions plotted in Figured] (open triangles) 
after removal of the estimated orbit. The errors bars shown in the bottom- left-hand corner 
represent the mean position error for the 35 epochs (see text), (b) As in (a) but now the sky 
distribution of the brightness-peak positions (plusses). To plot these positions, we added to 
the position in (a) for each epoch the difference between the position of the brightness peak 
(set i) and the position of the Gaussian or midpoint (set iii - see text). The mean position 
error in (b) is virtually identical to that shown in (a). The solid circle in each panel indicates 
the disk of the primary with a radius of 0.64 mas and centered at the origin by assuming that 
the center of the primary follows the estimated orbit. The dashed line passing through the 
center of each panel indicates the projected spin axis of the primary, assuming it is parallel 
to the orbit normal (see §EJ). The dashed circles in (b) separate the three regions used for 
the counting analysis described in § 14.4.31 
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using a bootstrap method (see lEfron fc Tibshiranilll993l ). More specifically, we regarded the 
35 residuals as the parent distribution and chose from this distribution a new set of 35 resid- 
uals, with each one being randomly selected from the parent distribution (with replacement), 
until we obtained 500 different sets of 35 residuals. For our estimate of the standard error 
of the rms scatter of the true distribution, i.e., before sampling by our observations, we take 
the standard deviation of the 500 values of rms scatter that we computed from the 500 sets. 
We used this bootstrap method to derive standard errors for all parameters estimated from, 
or compared to, the 35 residuals, both here and, where relevant, hereafter. 



4.4. Simulation of the Distribution of Position Residuals on the Sky and 

Comparison with Observations 

How can the scatter of the locations of brightness peaks be interpreted in terms of the 
geometry of the orbit and the relation to the primary? In other words, how closely is the 
seemingly preferred direction of the scatter of these locations related to the normal of the 
orbit and the spin axis of the primary, and how far from the surface of the primary do the 
brightness peaks originate? Motivated by our astrometric solution and the distribution of the 
brightness peak locations in Figure Ht, we constructed a three-dimensional model to simulate 
emission regions in the corona of the primary. We call this model the coronal emission model 
(CEM). We then considered different distributions of locations of emission regions to find 
the best match to our observations and thereby obtain a reasonable model for the location 
of the source of radio emission in IM Peg. 



4-4- 1- Latitude-independent Coronal Emission Model with Spillover Emission 

We start with the latitude independent version of the CEM. In this model we assume 
that on average the radio emission is centered, in projection on the sky, on the center of 
the stellar disk, and allow the location of peak brightness to fall with equal probability 
above any point on the stellar surface. We allow emission to occur at any height above the 
photosphere of the star, but assume that the probability that the emission peaks at a given 
height decreases exponentially with scale height, H® We then project the CEM onto the sky 



3 We are not asserting here that an exponentially decreasing emission probability is physically realistic. 
An exponential function provides via a single parameter an analytical means of estimating the statistical 
distribution of emission heights. Other functions (e.g., uniform emission probability) employ a sharp cut-off 
at an arbitrary height and, more importantly, were not able to reproduce in projection the full extent and 
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to compare our model distribution to the distribution of the locations of brightness peaks 
as plotted in Figure [3fr. For emission from locations not occulted by the disk of the star, 
we simply project that location onto the sky. For emission locations occulted by the disk of 
the star, we allow "spillover;" i.e., we move the predicted location of observed peak emission 
outward along the radial line on the sky to a location we can see slightly beyond the edge 
of the stellar disk (see below). Allowing for spillover is perhaps reasonable, since we had 
a 100% detection rate and observed extended radio structure at every epoch; however, as 
we shall later show, nature does not seem to have followed this path and there are other 
ways to understand the 100% detection rate. We do not account in any sophisticated way 
for the transmission of spillover emission through the corona, e.g., by considering scattering, 
absorption, and/or refraction along the line of sight. Instead, we arbitrarily allow a small 
radial distribution of this emission on the sky spanning a projected radius of 1.0-1.1 times 
the stellar radius. The adjustable parameter in the latitude-independent CEM is H. 

To execute the comparison, we computed several sets of 5000 locations each for the 
CEM, each set with the nominal value for the stellar radius but a different H. For each set 
we compared the rms of the distribution of the radio-emission locations to the 0.62±0.03 mas 
rms scatter of the observed residuals. We found agreement between the two rms values for 
H = 0.11 ± 0.03 mas; i.e., a scale height only 0.17 ± 0.05 times the 0.64 mas stellar radius. 
The standard error in H is the variation in that parameter needed to produce an increase 
or decrease in rms of the distribution of the radio-emission locations equal to the standard 
error (±0.03 mas) of the rms of the scatter in the residuals. We provide an illustration of 
the latitude-independent CEM in Figure Hfr. 

The figure illustrates that the latitude-independent CEM produces a circularly symmet- 
ric distribution of emission locations on the sky. Such a distribution fails qualitatively to 
describe the elongation in the scatter of the residuals. For a more direct comparison with the 
scatter of the 35 residuals, we produced several sets of 35 (random) emission locations from 
the latitude-independent CEM, and for none of them was there more than a 16% difference 
between the rms along the line at which it was greatest from that along the line for which it 
was least. In contrast, the residuals show a 52% ± 22% difference between the rms scatter 
along the line for which this rms scatter is maximum (p. a. = —38° ± 8°) and the line for 
which this scatter is minimum. 



distribution (in the three regions described in § I4.4.3[) of the observed residuals. 
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Fig. 4. — (a) Sky projection of 5000 random emission points from the latitude-independent CEM 
with spillover emission described in § 14.4.11 The scale height in the model is 0.17 times the 0.64 mas 
stellar radius, a value which produced the best match to the observed distribution. The dots 
represent emission points that are not occulted by the disk of the star. The "x's" (red in the 
colored version) just above the stellar disk represent spillover emission from points occulted by the 
disk of the star (see text), (b) As in (a) but now the points from the latitude-dependent CEM 
with spillover emission described in § 14.4.21 The scale height is 0.15 times the 0.64 mas stellar 
radius, again chosen so as to produce the best match to the observations, (c) As in (a) but now 
the points from the latitude-dependent CEM with reduced spillover emission described in § 14.4.31 
The scale height was 0.20 times the 0.64 mas stellar radius (again chosen for the best match), (d) 
Sky projection of 35 randomly selected emission points (plusses) from the latitude-dependent CEM 
with reduced spillover emission. For all panels, the solid circle gives the disk of the primary and 
the dashed circles separate the three regions used for the counting analysis described in § 14.4.31 



- 16 - 



4-4- 2. Latitude- dependent Coronal Emission Model with Spillover Emission 

To improve upon our CEM, we maintained our assumptions concerning spillover, but 
dropped our assumption of equal probability density per unit surface area. Instead, we al- 
lowed this probability density to vary with stellar latitude, A. This enhancement of our model 
was motivated by Doppler optical surface images of the IM Peg primary which show persis- 
tent (~l-3 years and possibly longer), hig h-intensity spot features i n both the mid-latitude 
range and directly oyer t he visible pole ( IBerdyugina et al.l l2000t iBerdyugina fc Marsden 



20061 ; iMarsden et al.l 120071 ) . Given the close alignment expected for the spin axis and the 
orbit normal (see §[3]), we assumed that the spin axis of the primary is inclined 73° to our 
line of sight, an angle equal to our VLBI-derived estimate of the inclination of the orbit (see 
Table |2]). We tried a number of different functional forms for the latitude dependence of the 
probability density distribution. We found that a distribution proportional to k + sin|A|, 
with k = 0.14^q;q4, yielded a good fit to the sky distribution of the residuals. For this k value, 
the mean probability density per unit surface area for emission near each pole (|A| > 70°) 
is 3.6^Qy higher than that for emission near the equator (|A| < 20°). The indicated ~70% 
confidence limits reflect the uncertainty of the elongation of the scatter of the residuals, as 
estimated with our bootstrap approach. 

We also varied the scale height, and found best agreement between the distribution of 
the CEM and the residuals for a value 0.15 ± 0.05 times the stellar radius. In Figure© we 
provide an illustration of the latitude-dependent CEM oriented on the sky with the spin axis 
lying along p. a. = —38°, to align with the axis of the elongation of the distribution of our 
35 residuals. 



4-4-3 ■ Latitude- dependent Coronal Emission Model with Reduced Spillover Emission 

The latitude-dependent CEM with spillover emission gives us a distribution of projected 
emission points which is consistent with the distribution of residuals both in overall extent 
and degree of elongation. To compare the distributions in more detail, we looked at the 
number of emission points and residuals in each of three regions at increasing radial distances 
from the center of the distribution. Based on the extent of the distribution, we chose the 
following regions: (1) <0.40 mas, (2) 0.40-0.80 mas, and (3) >0.80 mas. The regions are 
indicated by dashed circles in Figure [3^ and in each of the panels in Figure HI The results, 
given in Tabled show that the number of emission locations in region 2 is significantly higher 
for the latitude-dependent CEM with spillover emission than for the observed residuals. 
Figure shows clearly that, for the latitude-dependent CEM with spillover emission, many 
of the emission locations in region 2 are from spillover. Indeed, 35% of emission locations 
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in this region, and 28% of all points in this model are due to spillover. The percentage of 
instances for which the radio brightness peak arises from spillover emission may, however, 
be much lower. Here is why. In nine of 35 epochs, our images of IM Peg show two (or 
three) lo cal brightn ess maxima, with the maxima separated in four of those epochs by >0.5 



mas i see lPaper VIII ). Since simultaneous emission from multiple regions appears common for 
IM Peg, it is reasonable to expect that spillover emission be dominated in many instances 
by emission from the Earth-facing side of the primary. To investigate this possibility, we 
modified the selection of projected emission locations in the latitude dependent CEM with 
spillover emission. Instead of allowing each occulted location to become spillover emission, 
we allowed, first, only 50% of occulted locations to become spillover emission, replacing 
the other 50% with (randomly and independently drawn) non-occulted locations. We then 
decreased the fraction of allowed occulted locations to zero in intervals of 5% (i.e., 45%, 
40%, ...), adding the deleted ones at non-occulted locations. For each such fraction, we 
adjusted the scale height, H, so as to maximize the agreement between the distribution of 
the CEM emission locations and that of the residuals. We found the best agreement between 
the distribution of emission locations in the CEM and that of the residuals, considering 
the number counts in the three regions defined above, for a reduced frequency of spillover 
emission with a fraction of 10% ± 10% of allowed locations (see Table [3]), and a scale height 
0.20 ± 0.05 times the stellar radius. We provide in Figure Hb an illustration of this latitude- 
dependent CEM, with reduced spillover emission, oriented on the sky with the spin axis lying 
along p. a. = —38°. Note that for our best fit model, in which the spillover from 10% ± 10% 
of emission peaks in the occulted region yields an observed brightness peak, the percentage 
of all emission peaks which occur due to spillover is just 3% ± 3%. Thus, for a random 
sample of 35 emission locations, this value suggests that only <2 of the locations arise from 
spillover emission. We show one realization of 35 randomly chosen emission points from the 
latitude-dependent CEM with reduced spillover emission in Figure HJi. 



5. Discussion 



Our astrometric result shows that the sources of radio emission are consistent with 
their being centered on average on the primary, and also consistent with theoretical mod- 
els which propose that the radio emission in activ e close binaries is powered and confined 



by the magnetic field of the active star (see, e.g., iLestrade et al.l Il988l ; iMutel et al.l Il998 



Franciosini et al.lll999l ). Any constant or linearly changing offset of the center of the distri- 
bution of emission locations from the center of the disk of the primary over the 8.5 years of 
our VLBI observations is not necessarily expected to be zero. First, any latitude-dependent 
emission model when combined with shadowing of one of the pole regions due to the inclina- 
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Table 2. Orbit Parameters of IM Peg Radio Source 



Parameter Estimate Standard Error a 

— The parameters of the linear orbit model: — 
A sa (mas) -0.59 0.10 

A sS (mas) -0.66 0.11 

A ca (mas) 0.15 0.09 

A cS (mas) -0.23 0.11 

— The equivalent values of the more commonly used orbit parameters: — 
a (mas) 0.89 0.09 

Axial ratio b 0.30 0.13 
n{°) c 40.5 8.6 

T conj (JD) 2450342.56 0.44 

a See |Paper Vl for an explanation of our errors. 

b The axial ratio, i.e., the ratio of the minor to the major axis, is equal to the 
absolute value of the cosine of the inclination for our zero-eccentricity orbit. 
The axial ratio therefore corresponds to an inclination of 73 zfc 8° . 

c Position angle (east of north) of the ascending node. We follow the con- 
vention of identifying the ascending node as the one at which the source is 
receding from us as it passes through the plane of the sky. 



Table 3. Comparison of the Distribution of Residuals to the Distribution of CEM 

Emission Points 





Residuals 


lat.-dep. CEM with spillover a 


lat.-dep. CEM with reduced spillover* 3 




35 epochs 


35 points 


35 points 


Region 1 


6 ± 2 (17% ± 6%) 


2±1 (6% ±3%) 


5 ± 2 (14% ± 6%) 


Region 2 


22 ± 3 (63% ± 9%) 


29 ± 2 (83% ± 6%) 


24 ± 3 (69% ± 9%) 


Region 3 


7 ± 3 (20% ± 9%) 


4±2 (11% ±6%) 


7 ± 2 (20% ± 6%) 



a Latitude-dependent CEM with spillover emission as discussed in § 14.4.21 

b Latitude-dependent CEM with reduced spillover emission as discussed in § 14.4.31 The tabulated results 
correspond to a 10% allowed fraction of spillover points. 

°The value and standard error given for each region represent the mean number and standard deviation 
for 500 realizations of 35 randomly chosen CEM emission points (see text). 
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tion of the spin axis would predict a constant offset. Second, a systematic time- dependence 
of the latitude distribution of spot centers could arise given the apparent multi-year stellar- 



activity cycle of the IM Peg primary (e.g.. iBerdyugina et al.ll2000t IZellem et al.ll2010[ ). With 



the distribution of the radio emission locations at least partly linked (see below) to such spot 
features, an offset with a nonzero trend from 1997-2005 could plausibly contribute error to 
our proper-motion estimate. However, each of these pos sible caus es would not result in the 
offset being larger than the radius of the primary (see iPaper VI ). This estimate, together 
with the distribution's matching the disk of the primary, makes us believe that any offset at 
any time during our observations is indeed smaller than the radius of the primary. 

Our position residuals are scattered preferentially along an axis with p. a. of —38° ± 8°, 
which, to within the combined uncertainties, is equal to the p. a. of the sky-projected orbit 
normal of —49.5° ± 8.6°. Since the orbit normal is expected to be clos ely aligned with th e 
spin axis of the primary as for all synchronous RS C Vn systems (see IStawikowskil Il994l ) , 
the sources of radio emission appear to be linked to the spin axis of the primary. Indeed, 
a comparison of the scatter of the position residuals with our simulations shows that the 
probability density per unit surface area for radio brightness peaks is 3.6lo7 higher near the 
poles (A > 70°) than near the equator (A < 20°). Since Doppler images of IM Peg show 
the presence of persistent, high-intensity spot feat ures at the pole region of the primary 
(IBerdyugina &: Marsdenll2006l ; iMarsden et al.l 120071 ). our result provides statistical evidence 
that the radio emission regions occur at the same stellar latitudes as do active surface regions. 

Our simulations place restrictions on the altitude of the locations of radio-emitting 
structures in the corona of the IM Peg primary. In our best-fit model, ~2/3 of the brightness 
peaks arise within an altitude of just 0.25 times the stellar radius. Since the stellar magnetic 
field is presumably strongest near the surface, this result may not be surprising . However, it is 



differ ent from the pictures p resented for the close bi nary systems Algol (e.g., iLestrade et al. 



19881 ) and UX Arietis (e.g., iFranciosini et al.lll999l ). which show emission from magnetic- 
loop structures with heights greater than one stellar radius. If emission does occur high on 
magnetic loops for IM Peg, then most of these loops are small in height compared to the 
radius of the primary. 

Our simulations further allowed us to set at 6% the (la) limit on the likelihood (at any 
given epoch) that the brightness peak arises from spillover. Although this limit is model- 
dependent, we think that the more general conclusion, which is that the observed brightness 
peaks arise mostly from emission regions connected to the Earth-facing side of the primary, 
is robust. Spillover could, however, be responsible for lower-surface-brightness features of 
the radio structure. In many of our epochs, the IM Peg radio structure is multi-peaked or 
at least significantly elongated. While the highest peak in these instances is likely to be 
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associated with emis sion from t he Earth-facing side, spillover could very well contribute to 
the overall shape. In iPaper VIII . we present the full set of IM Peg images with the outlines 
of the disk and orbit of the primary superimposed, and discuss in detail the size and shape 
of the radio emission regions. 



6. Conclusions 

Here we summarize our results and give our conclusions: 

1. The sources of radio emission are on average located near the center of the disk of the 
primary, in particular being within 12 ± 15% of its radius of 0.64 ± 0.03 mas, provided any 
8.5-yr constant or linearly changing offset is (nearly) zero. There are theoretical arguments 
as well as observational evidence that any such offset is indeed smaller than the radius 
of the primary. Thus IM Peg is the second close binary with such (nearly) unambiguous 
identification. 

2. The positions of the sources of observed radio emission are scattered over an area on 
the sky slightly larger than the disk of the primary and preferentially along an axis with 
p. a. = —38° ± 8°. This axis is closely aligned with the sky-projected orbit normal (p. a. 
= —49.5° ±8.6°) and expected spin axis of the primary. 

3. Comparison of our observed positions with simulations suggests that the radio brightness 
peaks are more likely to occur at higher stellar latitudes than near the stellar equator, 
with the probability density per unit surface area being 3. 61° 7 times higher near the poles 
(A > 70°) than near the equator (A < 20°). The radio emission regions therefore show a 
dependence on stellar latitude similar to that exhibited by active regions on the primary's 
surface seen as dark spots with optical Doppler imaging. 

4. This comparison also suggests that these brightness peaks preferentially arise close to the 
surface of the primary, with ~2/3 of them located no more than 0.25 stellar radii above the 
surface. 

5. This comparison further suggests that the brightness peaks are mostly associated with 
emission regions on the Earth-facing side of the primary, with peaks caused by spillover from 
emission regions on the opposite side of the star occurring rarely, if ever. 
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